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Nearly twenty years ago, Beaurepaire and coworkers showed that when an uhrafast laser impinges on a 
ferromagnet, its spin moment undergoes a dramatic change, but how it works remains a mystery. While the 
current experiment is still unable to resolve the minute details of the spin change, crystal 
momentum-resolved techniques have long been used to analyze the charge dynamics in superconductors 
and strongly correlated materials. Here we extend it to probe spin moment change in the entire 
three-dimensional Brillouin zone for fee Co. Our results indeed show a strong spin activity along the A line, 
supporting a prior experimental finding. The spin active pockets coalesce into a series of spin surfaces that 
follow the Fermi surfaces. We predict two largest spin change pockets which have been elusive to 
experiments: one pocket is slightly below the A line and the other is along the A line and close to the L point. 
Our theory presents an opportunity for the time-, spin- and momentum-resolve photoemission technique. 

Crystal momentum-resolved techniques have been instrumental to understanding electronic and magnetic 
properties in a variety of materials from high Tc superconductors^ to magnetic materials. Angle-resolved 
photoemission (PE) is an example, where one can examine charge and lattice dynamics^ or pure electron 
dynamics in antiferromagnets^ along a few selected crystal momentum lines. This technique has now been 
extended to x-ray energy regions where its crystal-momentum resolving power is augmented by the element 
specificity. The first time-resolved photoemission measurement of the ultrafast spin dynamics in Co was demon- 
strated over a decade ago^, but it was not momentum-resolved. The situation did not improve much up to now^, 
with a few exceptions^. Majority of laser pulse durations are too long to fully resolve the spin polarization change 
at the earliest stage. Without any theoretical guidance, it is difficult to prove that the spin polarization really 
reflects the spin moment change in a sample. This is because experiments measure the spin polarization of those 
ejected electrons, not those electrons in the sample, while one is interested in the spin moment change of those 
electrons left behind. As a result, the response of spin-polarized photoelectrons to the laser pulse may differ from 
those electrons which contribute demagnetization in femtomagnetism. Weber et aV recently even showed 
different demagnetization times for the time- and spin- resolved photoemission and the time-resolved mag- 
neto-optics. This explains why Pickel et al.^ could only examine the momentum-resolved spin-orbit hybridiza- 
tion, not momentum-resolved spin moment changes which are a focus in femtomagnetism^"^°. 

Despite of those uncertainties in PE, crystal momentum-resolved techniques can provide a microscopic picture 
momentum by momentum and allow a direct comparison with theoretical calculations. Some of the latest 
experiments reported in a recent conference may represent a major breakthrough. Ross et aV^ showed that it 
is possible to enhance the experimental signal by an order of magnitude by collecting cascade electrons. In sharp 
contrast to almost all the prior studies, here the spin moment changes are crystal momentum-resolved, in a 
similar fashion as done for the charge dynamics. Undoubtedly, this will open a new frontier in femtomagnetism. 
A theoretical investigation is imperative. 

Face-centered cubic cobalt is a conventional 3 d ferromagnetic transition metal, and its magnetic properties 
have been thoroughly studied over two decades^^. In the realistic calculation, we use the experimental lattice 
constant of 3.549 A^^ The obtained magnetic spin moment of Co is about 1.66 fi^, which agrees well with the 
previous theoretical value 1 .64 ^b^^- While more advanced theoretical methods are available for smaller systems^^, 
at present it is difficult, if not impossible, to use them to compute dynamical processes for a solid. Our theoretical 
calculation is based on the density functional theory (DFT) using the full potential augmented planewave method 
as implemented in the WIEN2k code^^'^^. For the exchange-correlation functional we use the generalized gradient 
approximation (GGA) in the form of Perdew-Burke-Ernzerhof (PBE) functionaP^. Our spin moment and band 
structure fully agree with the previous experimental results and theoretical results^^'^^"^\ In the laser excitation 
both the Kohn-Sham (KS) ground and excited states are important; to accurately take into account the delocalized 
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nature of excited states, a Muffin-tin radius R^t is set to 2.36 a.u. and 
a large product of the Muffin-tin radius and the planewave cutoff, i.e., 
RmtKmax^ is Set to 9.5. Thus, more basis functions are included to 
accurately describe high lying states at the laser photon energy. The 
number of the basis functions is 95. After solving the KS equation 
self- consistently, we compute the optical transition^^'^^ and spin mat- 
rices'^, and get ready for the laser-induced spin dynamics. 

Although significant work has been devoted to the understanding 
of ultrafast demagnetization, for example electron-phonon scatter- 
ing theory^^'^^, superdiffusive spin transport mechanism^^'^^, and 
others^^~^°, the microscopic origin is still unclear. Alternatively, we 
propose a demagnetization mechanism based on the electronic tran- 
sitions, where the spin- orbit coupling (SOC) plays a role in the spin 
moment change. Under the influence of a laser pulse, electrons are 
excited out of the Fermi sea. The entire electron dynamics can be 
quantitatively computed through the Liouville equation for the 
density matrix p at each ih(dpj^/dt)=[Ho-\-Hi,pj^]. This is a 
coupled equation involving transitions between all the band states 
and is solved numerically. Here Hq is the system Hamiltonian, and 
Hi is the interaction between the laser and the system, i.e., 
Hi= ^^^^ ,.Dk-ij-E{t)cl.Ck,jy where D is the dipole operator, and 

is the electron creation (annihilation) operator for band / 
at k. The laser field is E{t) = eAo exp [ — f / t^] cos{cDt), where e is the 
laser polarization direction, Aq is the field amplitude, co is the laser 
frequency, t is the time and t is the laser pulse duration. Once we 
obtain the density matrix, the time- dependent spin moment can be 
computed by tracing over the product of the density matrix and spin 
matrix. The major strength of our simulation is that we realistically 
include the interaction between the laser system and the ferromag- 
netic sample. There is no fitting parameter in our model. The major 
weakness is that the electron- electron interaction is not included. 
There are two significant findings of our study. First, although there 
are many excitations among electrons, only those electrons, whose 
initial and final spin moments are different, can contribute strongly 
to the total spin moment change. This process is mediated by SOC. 
Second, we demonstrate that there are large pockets of the isosur- 
faces of spin moment change in the Brillouin zone, which are access- 
ible to the time- and crystal momentum- resolved photoemission 
experiments. 

We emphasize that the dipole matrix elements are not updated 
during our time dependent simulations. This is different from the 
standard TD-DFT where the dipole element is updated dynamically. 
Unfortunately, a TD-DFT calculation is extremely difficult if it is not 
impossible because nearly half million k points are needed to con- 
verge the spin moment change'^. This situation is unlikely to change 
in the near future. Our method relies on the ground state property of 
DFT, and is a compromise between the accuracy and the efficiency. 
Naturally, we are fully aware of the limitation of our formalism. First, 
its validity has been checked against a correlated modeP^. Second, we 
also checked against the quantum chemistry calculation in a small 
system'^, which is carried out by the Gaussian package. The agree- 
ment is very good. Third, we carried out a calculation using the TD- 
DFT in a small system and then compared our results. We find that 
the major difference in the electric polarization appears at the earliest 
time, where the exchange kernel undergoes the major change. The 
difference is that in our method there is a smaller delay between the 
laser peak and the system's response, but the TD-DFT has a slightly 
larger delay (larger than 1-5 fs). All the subsequent evolution is 
nearly identical. Since the spin change occurs on a few hundred fs, 
this short delay does not affect our results significantly. When a laser 
pulse becomes too strong, our formaUsm becomes less accurate, since 
the laser perturbs the system strongly. 

It is well known that photon energy plays a vital role in femto- 
magnetism. This is because it acts as a selector of electron transition 
window, which subsequently dominates the spin moment change 



associated with bands. Thus, we optimize the photon energy where 
the pulse duration of 12 fs and the field amplitude of 0.05 V/A are 
fixed during the optimization process. Our results show that the spin 
moment change depends critically on the photon energy. When the 
photon energy is about 1.9 eV, the laser pulse generates the max- 
imum spin moment reduction. The obtained demagnetization curve 
is similar to that obtained in fee Ni^^ (not shown for brevity). The spin 
moment first drops sharply and then reaches its minimum at 23.6 fs. 
After that the spin recovers itself and oscillates around an average 
value of 1.638 ^b- In comparison with the result in fee Ni, we notice 
that the first minimum appears at a later time which originates from 
the weaker SOC in fee Co. 

In order to reveal further insights into the spin moment change, we 
disperse the time-averaged spin change (after 90 fs) in the first 
Brillouin zone. While there are a few isolated hot spin spots^^, mojor- 
ity of the spots with large spin moment change coalesce into larger 
segments. Their intensity map is defined as the spin surface, in ana- 
logy to the Fermi surface. There are three reasons why these spin 
surfaces exist. First, due to the laser photon energy, only a limited 
number of unoccupied conduction bands are accessible from the 
occupied states, the top of which are at the Fermi level. In our case, 
these bands appear about 2 eV (our photon energy) above the occu- 
pied band states. Those bands form a surface. Second, among those 
unoccupied bands, a small portion is excited strongly because their 
transition energy matches the photon energy. Third, among those 
excited bands, only those bands that have a large spin moment 
change can show up in the intensity map or the spin surface. Since 
experiments are almost always performed along some high- sym- 
metry lines and planes, we slide through several traditional high 
symmetry surfaces. Figure 1 displays the spin surfaces within four 
planes, where the false colors represent the amount of spin moment 
changes with the color bars depicting their true values. In Figs. 1(a)- 
1(d), five curves numbered from 1 to 5 represent their Fermi surfaces. 
Our results are very insightful. There are five spin surfaces that have 
an appreciable spin moment change; and other regions have a much 
smaller spin change. We start with the F-K-U-W-K plane (see 
Fig. 1(a)). Within this plane, the spin surface closely follows the 
Fermi surface (see Fermi surface 1 in the figure), but there is a 
distance between these two surfaces. This distance is a result of the 
laser photon energy Hco. Depending on Hoj, the spin surface can move 
away from or toward the Fermi surface. From the color, one can tell 
that the upper part of the spin surface is spin-moment-enhanced, 
while the lower part is spin-reduced. Within the same plane, another 
spin surface follows Fermi surface 2. This surface has the largest spin 
reduction in the entire BZ. The third spin surface follows Fermi 
surface 3, but the amount of the spin reduction is smaller than the 
first two. 

Experimentally, Pickel et al.^ did detect the spin activities along the 
A direction. This is the predictive power of our theoretical method. 
However, Pickel et could not attribute those spin activities 
directly to the demagnetization, due to the reasons as discussed in 
the Introduction. Nevertheless, our theory does corroborate the spin 
change along the A direction. More importantly, our result reveals 
two much larger spin change pockets: one is slightly below the A line 
and close to the X point and the other is the largest spin change along 
the A direction and closer to the L point. Therefore, experimentally 
there is an opportunity to detect stronger spin signals. Future experi- 
ments can directly test our prediction. 

Pickel et aV also theoretically showed additional spin hot spots 
along the E line, but this is not observed in our simulation. Instead, 
our data shows that there is no major spin change along that dir- 
ection. This demonstrates that while a strong spin-orbit hybridiza- 
tion is important to the spin response, it is not a sufficient condition 
of a spin moment change. We also notice that in comparison with our 
band structure (see the bottom panel of Fig. 1) and prior experi- 
mental'^ and theoretical results'^'^°, their band structure is unexpect- 
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Figure 1 | (a) Spin moment change AmJ in the T-X-W-K and T-X-U-L-K 
planes. The intensity map denotes the spin moment change (in units of 
^b)- The curves denote the Fermi surfaces (same for all the other figures). 
The spin surfaces follows the Fermi surfaces, (b) Spin moment changes in 
the L-K-W-U-W plane are concentrated around Fermi surface 4. (c) Spin 
reduction in the X-W-W plane is the smallest among those four planes, 
(d) Spin moment change in the X'-W-W plan. The spin surface is away 
from the Fermi surface. Bottom panel: Band structure of fee Co. 



tions, 



^ijkl- 



In order to understand the above results, we first note that since the 
photon energy is fked, there are only a limited number of the states that 
can be excited strongly. This can be seen from the band structure itself 
(bottom panel of Fig. 1). For instance, close to T point, transitions are 
difficult since there is no state available around 1.9 eV. Second, without 
SOC, the optical transition conserves the spin and no spin change should 
be observed. This also applies to the Coulomb and exchange interac- 

where is the distance between two electrons. Electrons can stay 
in any orbitals / through /, but they can only carry two different spin 
indices Gi and (72. The fundamental reason is because the electro- 
static interaction does not operate on the electron spins. 
Mathematically, V^e commutes with the total spin^^. This has an 
important consequence on the Stoner excitation. Although the 
Stoner excitation is frequently cited as a way to the demagnetization, 
its demagnetization also relies on a spin symmetry breaking term, 
which is often SOC. However, the traditional theory hides this origin 
by invoking a thermal reservoir^. Third, with SOC, spin moment 



(a) 




(b) 



f 



edly upward-shifted by about 0.3 eV or the Fermi level is down- 
shifted by the same amount. At present, we do not know the reason 
behind this discrepancy. 

Figure 1(b) shows the spin dispersion in the L-W'-U-K plane. Due 
to SOC, W and W points become nonequivalent. The major spin 
changes appear in two different regions. The first one is around the L 
point, and the second one is around the W point. Different from 
Fig. 1 (a), here the spin reduction and spin enhancement appear in the 
same region. However, the spin reduction surface is closer to Fermi 
surface 3, but the spin enhancement surface is further away from 
Fermi surface 3. On the other hand, the spin surface around the W 
points only shows the spin decrease. In general, this plane has more 
points with spin moment reduced than those with spin enhanced. 

Figure 1(c) illustrates the spin change in a plane that is perpen- 
dicular to the A line. The spin surface (a square) differs from the 
Fermi surface 5 (a circle). The spin change is the smallest among 
these four planes. Different from Fig. 1(c), Fig. 1(d) shows a stronger 
spin reduction. This plane is perpendicular to the axis. The main 
changes are closer to the W points. 



(c) 




Figure 2 | (a) Average spin moment change (averaged after 90 fs) is 
dispersed in the Brillouin zone. The isosurface value is taken at —0.3 jib- 
The high-symmetry points and lines are labeled, (b) and (c) Two selected 
Fermi surfaces which are related to the isospin surfaces closely. 
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changes become possible. Our data shows that the most important 
contribution to the spin moment change comes from those transi- 
tions with a bigger difference between initial and final spin moments. 
The challenge is that these transitions are not optically favorable. In 
other words, transitions with a large spin change are often weakly 
dipole- allowed. On the other hand, those strongly dipole-allowed 
transitions have smaller spin changes. As a result, the true demag- 
netization as seen in these planes is a compromise between these 
two competing effects, which is characterized by the spin- dipole 
factor^^'^^ 

Figure 2(a) shows the three-dimensional spin moment dispersion 
in the crystal momentum space. This is an isosurface with spin 
moment change of —0.3 which contributes demagnetization. 
This dispersion is constructed from 73763 irreducible k points, 
among which only 681 k points with the spin moment change are 
smaller than —0.15 fi^ (more negative). For our present laser para- 
meters, the largest moment change is —0.72 at (46, 44, —44)/ 104 
(in the unit of 2n/a, where a is the lattice constant of Co). If all the 
other k points had the similar spin reduction, this would mean 56% 
reduction in the spin moment. Unfortunately, a large portion of k 
points are silent. In addition, not all the k points have their spin 
moments reduced, and some k points see an increase in spin 
moment. For instance, the largest increase in spin moment is 
+ 0.44 ^B at (76, 22, 0)/104. This decreases the amount of the spin 
reduction. The key message from this figure is that the spin change is 
not evenly distributed and is well localized behind the L-U-W and the 
hexagonal planes which are consistent with our above discussion. To 
help locate these spin isosurfaces, in Figs. 2(b) and 2(c) we show the 
corresponding Fermi surfaces. They closely resemble each other. 
Future experiments can use the Fermi surface as a guide to detect 
those spin isosurfaces. If these isosurfaces are measured as a function 
of time, one can construct a dynamic three dimensional crystal- 
momentum-resolved dispersion of the spin moment change, a movie 
in real time domain. Our Supplementary Material shows our theor- 
etical movie. An experimental verification is much needed. 

In conclusion, for the first time the ultrafast spin moment change 
is dispersed in the entire Brillouin zone for fee Co. Our theory finally 
links those spin activities observed experimentally to the genuine 
spin moment change. Our results support the experimental obser- 
vation of the strong spin activities along the A line^. More impor- 
tantly, we reveal two even stronger but previously unobserved spin 
hot pockets. One pocket is slightly below the A line and the other is 
along the A line. This presents an excellent opportunity for future 
experiments. Further developments in experimental techniques will 
finally allow one to probe the dispersion of spin moment dispersion k 
point by k point, a necessary step to fully understand femtosecond 
magnetism. 
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